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We present a measurement of the tt pair production cross section in pp collisions at yfs = 1.96 
TeV utilizing approximately 425 pb~ of data collected with the DO detector. We consider decay 
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channels containing two high pr charged leptons (either e or /i) from leptonic decays of both top- 
daughter W bosons. These were gathered using four sets of selection criteria, three of which required 
that a pair of fully identified leptons {i.e., e/j., ee, or fi/j,) be found. The fourth approach imposed 
less restrictive criteria on one of the lepton candidates and required that at least one hadronic jet 
in each event be tagged as containing a b quark. For a top quark mass of 175 GeV, the measured 
cross section is 7.4 ±1.4 (stat) ±1.0 (syst) pb and for the current Tevatron average top quark mass 
of 170.9 GeV, the resulting value of the cross section is 7.8 ±1.8 (stat+syst) pb. 

PACS numbers: 13.85.Lg, 13.85.Qk, 14.65.Ha 



I. INTRODUCTION 
A. The Top Quark 

The top quark, first observed by the CDF and DO col- 
laborations in 1995 P, is the heaviest elementary par- 
ticle so far observed. Its mass is sufficient to allow decay 
to hypothesized particles such as the charged Higgs and 
to probe electroweak symmetry breaking physics. At the 
Fermilab Tevatron Collider, top quark production occurs 
predominantly in top-antitop quark (tt) pairs. For a cen- 
ter of mass energy of y/s = 1.96 TeV, leading order QCD 
suggests that tt production results from quark-antiquark 
annihilation about 85% of the time, while gluon-gluon 
fusion is responsible for the remaining 15% [3j]. Recent 
theoretical calculations predict, for an assumed top quark 
mass (ffl ( ) of 175 GeV, an inclusive top quark pair pro- 
duction cross section at y/s = 1.96 TeV of 6.7 pb with 
an uncertainty of less than 15% 0, If the observed 
production cross section were to differ significantly from 
the standard model prediction, it would be evidence of 
new physics, such as exotic top quark decays or new pro- 
duction mechanisms such as tt resonances [6j . Significant 
deviation among measured cross sections obtained from 
the observations of different top quark decay channels 
would also indicate the presence of new physics. It is 
therefore important to precisely measure the top quark 
pair production cross section using each possible final 
state. Previous measurements by the CDF and DO ex- 
periments [1, 0, H, Q show good agreement with the the- 
oretical expectation within uncertainties. The most pre- 
cise cross section measurement reported by DO is 6.6±1.0 
pb @ in the lepton±jets final state and using secondary 
vertex tagging algorithm to identify b jets. 



B. Top Quark Decays and the Dilepton Signature 

According to the standard model, the top quark decays 
almost 100% of the time to a W boson and a b quark. 
For approximately 6% of tt pairs, both W bosons decay 
leptonically to generate a final state containing a pair of 
electrons, a pair of muons, or an electron and a muon ||. 
This produces a unique event signature consisting of two 
high transverse momentum (j>t) charged leptons, signifi- 
cant missing transverse energy ($t) from the associated 
neutrinos, and two high pt jets from the b quarks. 

Despite low branching ratios relative to channels with 



hadronic W boson decays, the dilepton channels are ad- 
vantageous for study because few standard model back- 
ground processes have two high px leptons and neutrinos 
in their final states. Those which do usually do not con- 
tain two high pt jets. For example, electroweak diboson 
production can result in two isolated, high px leptons 
and neutrinos, but suffers from a low cross section and 
can be discriminated against by requiring high pt jets. 
Drell-Yan production of (Z/7*)±jets events has no direct 
decay process to dilepton final states with real neutrinos. 
(Z/j*) decay to r particles produces neutrinos but suf- 
fers from a low branching ratio and a softer lepton px 
spectrum relative to top quark events. 



C. The Content of this Article 

This paper describes a new measurement of top quark 
decays to final states containing a pair of electrons or 
muons, or one electron and one muon. Section [TTJ con- 
tains a description of the experimental setup used to col- 
lect the data used for the measurement. A discussion of 
the Monte Carlo samples that aided our interpretation of 
this data is in Sec. IIIII Section ITVl includes a description 
of the triggering system used to acquire the data, and 
Sec. IVl contains descriptions of the offline reconstruction 
techniques used to compute the physical quantities criti- 
cal to the extraction of the top quark signal. Discussion 
of the methods used to identify each dilepton decay mode 
in the data sample is in Sec. lVIl Finally, the computation 
of the top quark pair production cross section is described 
in Sec. IVIII and the result is summarized in Sec. IVIIII 



II. EXPERIMENTAL APPARATUS 

A. The Fermilab Tevatron Collider's Run II 

The Fermilab Tevatron Collider, a proton anti-proton 
accelerator, collided beams at a center-of-mass energy of 
1.8 TeV during the period of operation (Run I) between 
1992 and 1996. The DO detector, one of two multipurpose 
detectors designed to study the high energy collisions at 
Fermilab, collected approximately 120 pb -1 of data dur- 
ing Run I [10[ . After significant improvements to both 
the accelerator and the DO detector, Run II began in 
March 2001 with the collider operating at a center-of- 
mass energy of 1.96 TeV. The increased energy brought 
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an increase in the top quark pair production cross section 
of «30%. The analyses discussed in this paper are based 
on approximately 425 pb -1 of data collected by DO be- 
tween April 2002 and August 2004. DO has performed a 
similar measurement using ^230 pb" 1 of data [8|. 



B. The DO Detector 

This section presents an overview of the experimental 
apparatus, emphasizing the subsystems most relevant to 
the ti production cross section measurement. A more 
complete description of the upgraded experiment can be 
found in Ref. [111 ]. 

The DO detector comprises three major subsystems 
which together identify and measure the energy or mo- 
mentum of electrons, jets, muons, and (indirectly) neu- 
trinos — all of which can be found in the final states of 
ti decays. The subsystems are the central tracking de- 
tectors, a uranium/liquid-argon calorimeter, and a muon 
spectrometer. 

The spatial coordinates of the DO detector are defined 
as follows: the positive z direction is along the direction 
of the proton motion while positive y is defined as up- 
ward with respect to the detector's center, which serves 
as the origin. The polar angle 9 is measured with respect 
to the positive z direction and the azimuthal angle <fi is 
measured with respect to the positive x direction. The 
radial distance r is the perpendicular displacement from 
the z axis. The polar direction is more usually described 
by the pseudorapidity, defined as 77 = — ln(tan0/2). 

The central tracking detectors consist of a silicon mi- 
crostrip tracker (SMT) and a central scintillating fiber 
tracker (CFT) located within a 2 T solenoidal magnetic 
field. Together these detectors are responsible for locat- 
ing the position of the hard scatter and for measuring the 
trajectories and momenta of charged particles. The SMT 
can also locate displaced, secondary vertices which aid in 
heavy quark tagging. It is composed of high-resistivity 
silicon sensors arranged in barrels and disks to maximize 
the detector surface area perpendicular to charged parti- 
cle trajectories. The barrel detectors provide tracking in- 
formation at central values of (< 1.5), while the disks 
extend coverage out to \rj\ « 3.0. The CFT is constructed 
from scintillating fibers mounted on eight concentric sup- 
port cylinders. Each cylinder supports an axial layer of 
fibers oriented along z and a stereo layer oriented at a 
slight angle with respect to z. The outermost cylinder 
provides coverage for |?y| < 1.7. 

The liquid-argon calorimeter surrounds the central 
tracking detectors. In addition to providing energy mea- 
surements for electrons, photons, and jets, it can dis- 
tinguish showers generated by electrons or photons from 
those produced by hadrons. The calorimeter also plays a 
critical role in the measurement of the event-wide trans- 
verse energy balance used to identify neutrinos. The sys- 
tem is composed of three parts: a central calorimeter 
(CC) which provides coverage to \rj\ « 1 and north and 



south endcap calorimeters (EC) which extend coverage 
to \tj\ ~ 4. Because each calorimeter is housed in its 
own cryostat, there is a gap in coverage between each 
EC and the CC, the region defined by 1.0 < \rj\ < 1.4. 
To partially compensate for this, an intcrcryostat detec- 
tor (ICD) made of a series of scintillating tiles is located 
between the CC and EC cryostats. 

Each calorimeter section has three subsections: an in- 
ner electromagnetic (EM) section which uses thin ura- 
nium absorber plates, a fine hadronic section which uses 
uranium-niobium alloy plates, and a coarse hadronic sec- 
tion which uses copper or stainless steel absorber plates in 
the CC or EC, respectively. The calorimeters are trans- 
versely divided into projective towers, so-called because 
the rays along which the calorimeter cells are oriented 
project outward from the interaction center. Each tower 
layer is further divided into segments of size Arj x A<^> = 
0.1 x 2-7r/64, except for the third layer of the EM sec- 
tion which is segmented twice as finely to allow for more 
precise measurement of the EM shower centroid. 

The muon spectrometer surrounds the calorimeter 
cryostats and uses a combination of wire chambers and 
scintillation counters to obtain precise muon spatial and 
timing information, respectively. Like the calorimeter, 
the muon spectrometer consists of three separate sub- 
detectors. The central detector covers approximately 
1 77 1 < 1.0 and the forward systems extend to \r/\ 2. Each 
system contains three layers of instrumentation, and a 1.8 
T iron toroidal magnet is located between the innermost 
and second layers. Each layer contains both wire cham- 
bers and scintillation counters. The scintillators have 
response times sufficiently fast to allow for both muon 
triggering and out-of-time background rejection. 

The wire chambers in the central region are propor- 
tional drift tubes (PDTs) oriented to provide maximum 
resolution for measuring muon bending angles produced 
by the toroidal magnetic field. The innermost central 
scintillation counters are segmented in 4.5° increments in 
(f> to match the CFT segmentation. Each layer of the for- 
ward spectrometers contains several strata of mini drift 
tubes (MDTs) and a set of scintillation counters referred 
to as pixels. The pixels are projectively arranged from 
the interaction point with a <f> segmentation of 4.5° and 
an 77 segmentation of «0.12. 

The luminosity measurement is based on the rate of 
inelastic pp collisions observed by the luminosity mon- 
itors (LM) mounted in front of the EC cryostats at 
z = ±140 cm. The LM consists of two arrays of 24 
plastic scintillator counters with photomultiplier read- 
out, and covers the \rj\ range between 2.7 and 4.4. The 
uncertainty on the luminosity measurement is currently 
estimated to be ± 6.1% [13]. 



III. EVENT SIMULATION 

Selection efficiencies for ti signal events and back- 
ground survival rates for each of the analyses were com- 
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puted using Monte Carlo simulations of each of the 
physics processes contributing to the observed event 
yields. This section provides some details regarding the 
generation of the Monte Carlo samples used. 

Simulation began with initial parton generation. In 
general, this was achieved using the ALPGEN [l3[ genera- 
tor, which contains the exact leading-order (LO) matrix 
elements for the processes discussed in the following sec- 
tions. Unless otherwise specified, output from ALPGEN 
was then convoluted with the CTEQ5L [3| parton dis- 
tribution functions (PDFs). Parton showering was car- 
ried out using pythia flal- Decays of B mesons were 
simulated with evtgen llal and r-lepton decays were 
simulated with tauola [17j . 

After the modeling of quark and gluon hadronization 
and unstable particle decays, the list of gene rated ob- 
jects was passed through a GEANT-based [18| model of 
the DO detector. This provides a detailed simulation of 
the effects of detector composition and geometry. Res- 
olutions for momenta and energies of leptons and jets, 
as well as efficiencies for their identification, were de- 
termined in data and compared to their counterparts in 
Monte Carlo. Observed discrepancies were used to cor- 
rect the simulated samples. 

A. tt Production 

Detector acceptance, object reconstruction efficiencies, 
and the effects of kinematical cuts were estimated with a 
sample of simulated tt — > 11 + X decays, where I = e, (jl, 
or t. Seven samples were generated with the following 
values of top quark mass (nit): 140, 160, 175, 190, and 
210 GeV. These were used to parameterize the signal 
acceptances as functions of m t . The central value of the 
cross section was computed for m t — 175 GeV. 

B. (Z/7*)+jets Processes 

The largest background to the tt dilepton signal arises 
from Drell-Yan Z/j* production and leptonic decay with 
associated production of one or more jets. To aid in the 
estimation of these backgrounds, we generated Zj^f* — > 
U events with one or two partons. For each lepton fla- 
vor we generated three M« regions: 15-60 GeV, 60-130 
GeV, and >130 GeV. The relative weights of the three 
samples were determined from the ratios of their LO cross 
sections. 

The absolute normalizations of these background sam- 
ples were set by the number of Z/7*+jets — > £€+jets 
events observed in control samples selected from data. 
These were chosen by requiring that the reconstructed 
dilepton mass be near the Z boson mass so that the 
samples were rejected by the signal selection criteria de- 
scribed in Sees. IVI CI and IVI Dl This normalization was 
carried out at an early stage of event selection where the 
tt yield is a negligible fraction of the selected sample. 



The efficiency of further kinematical selections for Z/^j* 
events was then derived from the details of the simulated 
samples. 

C. Diboson+jets Production 

WW -> 11 + X and WZ -> 11 + X (where I = e or p) 
production in association with one or two jets contributes 
to a lesser degree to the selected samples. As for all 
the other Monte Carlo samples used, we generated these 
processes at LO with ALPGEN, but the impact of PDFs 
was simulated using CTEQ4L[19|. The resulting samples 
were normalized using the ratio of NLO to LO diboson 
production cross sections calculated without explicit jet 
requirements. (2p| . 

IV. TRIGGERING 

The analyses described in this paper made use of data 
collected by triggering on the presence of objects consis- 
tent with the dilepton signature: electrons, muons, cen- 
tral tracks, and jets. Correlations between these objects 
and event-wide variables like $t, though available at all 
trigger levels, were not utilized in data collection. This 
section begin with a brief description of the DO trigger 
system and then provides a description of the triggering 
conditions used to collect the dilepton samples analyzed. 
A more detailed discussion of the triggering system is 
available in Ref. [li"l |. 

A. The DO Triggering System 

The DO triggering system consists of three separate 
levels, each of which examines successively fewer events in 
ever greater detail. The first stage (Level 1) is a collection 
of custom hardware triggers that accepts data from all 
the major detector subsystems at a rate of 1.7 MHz and 
generates an acceptance rate of around or below 2 kHz. 
In the second stage (Level 2), microprocessors associated 
with each detector subsystem reconstruct physics objects 
which are passed on to a global processor that generates 
decisions based upon all the objects in an event. Level 2 
provides a maximum accept rate of around 1 kHz. The 
final trigger stage (Level 3), applies more sophisticated 
algorithms to data from precision readout of the detector 
components to further reduce the overall acceptance rate 
to around 50 Hz. Events passing Level 3 are written to 
tape. 

At Level 1, the muon trigger searches for patterns of 
scintillator and wire chamber hits consistent with muons 
traversing the multiple layers of the muon detector. Loose 
Level 1 muons are constructed from scintillator hits only, 
while tight muons include corresponding patterns of hits 
in wire chambers [29j |. Additionally, some Level 1 muon 
triggers require that a matching track be found by the 
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central track trigger (CTT). CTT tracks are found by 
analyzing patterns of axial CFT hits. All eight axial lay- 
ers must register a hit, and the curvature of the resulting 
patterns provides a pt estimate that is used for a thresh- 
old requirement. 

At Level 2, the muon- finding algorithm uses more pre- 
cise timing information to improve the quality of muon 
candidates. In general, a combination of wire and scin- 
tillator hits both inside and outside the toroid iron is 
required. Level 3 uses tracks found in the central tracker 
to identify the most probable position for the hard scat- 
ter. This position, also called the primary vertex, is used 
to refine momentum estimates from reconstructed muon- 
track bending by the toroidal field, and the result is used 
to apply momentum threshold requirements. Addition- 
ally, Level 3 is capable of reconstructing central tracks 
with hits missing and its algorithms make use of CFT 
stereo information. 

The Level 1 calorimeter trigger inputs are electro- 
magnetic (EM) and hadronic (H) trigger tower energies 
summed over a transverse area of Arj x A<fi = 0.2 x 0.2. 
Electron candidates only include energy collected in the 
EM section of the calorimeter, while jet candidate ener- 
gies include the H towers. At Level 2, calorimeter objects 
are reconstructed from trigger towers using the Level 1 
objects as seeds. The Level 2 jet algorithm clusters 5 x 
5 groups of towers centered on the seed towers. Electron 
candidates are formed by clustering each EM seed with 
the highest Et neighboring tower. 

The Level 3 calorimeter triggers use the precision read- 
out chain and the reconstructed primary vertex posi- 
tion to improve energy and position resolution relative 
to Level 2. Jets and electrons are formed using a simple 
cone algorithm 21]. Loose Level 3 electrons must have 
most of their energy deposited in the EM layers and they 
must meet basic shower shape criteria. Tight Level 3 elec- 
trons must survive additional shape criteria. Additional 
background suppression is also achieved in some triggers 
by requiring that a matching central track be found. 



B. Dilepton Triggers 

The triggers used to collect the dilepton samples re- 
quired that the lepton signatures distinguishing each 
channel were present at multiple triggering levels. In or- 
der to reduce rate but maintain overall trigger efficiency 
for a given channel, a logical OR of multiple triggers 
having different conditions tightened in a complemen- 
tary manner was sometimes used. Brief summaries of 
the trigger conditions used for each analysis channel are 
presented here, and more detailed breakdowns of the re- 
quirements are available in Appendix |A"1 

The e/i triggers required that an electron with an Et 
of at least 5 GeV and a loose muon were found at Level 
1. In some cases a Level 2 muon was also required, but 
otherwise the remaining conditions involved electrons re- 
constructed at Level 3. A loose Level 3 electron with 



Et > 10 GeV was always included in the requirements, 
and at higher luminosity the energy threshold was in- 
creased and this requirement was combined in an OR 
with a tight electron with Et > 5 GeV. 

The dielectron triggers usually included the require- 
ment that two electrons, each with Et > 6 GeV, were 
found at Level 1. In some cases only one electron with 
Et > 11 GeV was required at Level 1. A Level 2 re- 
quirement was only included for later periods containing 
high luminosity conditions. It required that two elec- 
trons, each with Et > 18 GeV, be found. The Level 3 
condition always included at least one loose electron with 
Et > 10 GeV. For later data taking periods, a second 
electron was added to the Level 3 condition, and the en- 
ergy threshold and quality requirements were tightened. 

To maximize efficiency, the /i/i channel made use of 
high pt single muon triggers and switched to dimuon 
trigger requirements when the single muon triggers were 
prescaled due to high rates. The single muon triggers 
required a tight muon at Level 1, while the dimuon trig- 
gers used loose Level 1 muons. All triggers used for the 
dimuon channel demanded that one muon be found at 
Level 2, sometimes with a pt > 3 GeV requirement. 
The Level 3 requirements also involved single muon sig- 
natures, but complementary conditions were sometimes 
combined in a logical OR. These signatures included a 
Level 3 muon with a pt of at least 6 GeV or a Level 3 
central track with a pt of at least 5 GeV. 

Since the ^+track channels did not require that two 
identified leptons be found in each candidate event, they 
relied on high pt single lepton triggers. In some cases 
these triggers included jet requirements. At Level 1, the 
e+track triggers demanded the presence of either one EM 
object with Et > 10 GeV or two EM objects, each with 
Et > 3 GeV. In some cases the single electron condition 
was coupled with the requirement that two Level 1 jets 
were found, each with Et > 5 GeV. The Level 2 condi- 
tions included the presence of one electron with Et > 10 
GeV. Some triggers also asked that two jets, each with 
Et > 10 GeV, be found at Level 2. Level 3 requirements 
included at least one electron with Et > 15 GeV. Some 
triggers also required that two jets, each with Et > 15 
GeV, be found at Level 3. 

The /x+track triggers required that at least one loose 
Level 1 muon, sometimes with a matching central track, 
be found. Some triggers also demanded that at least one 
jet with Et > 3 GeV be found. The Level 2 conditions 
usually required that one muon be found and sometimes 
also required one jet with Et > 8 GeV to be present. 
Level 3 conditions alternately included 1 jet with Et > 
10 GeV, one muon with pt > 15 GeV, or one central 
track with pt > 10 GeV. Sometimes the track and muon 
requirements were combined in a logical OR. 

The efficiencies of trigger conditions on single objects 
were estimated using data samples selected to remove 
triggering bias. The efficiency for a tt event to satisfy 
a trigger condition was then estimated by folding per- 
muon, per-electron, and per-jet efficiencies into Monte 
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Carlo simulated events. A similar process was used for 
those background estimations that are based upon simu- 
lation. 

Trigger terms related to electrons, muons, and tracks 
were analyzed using reconstructed Z boson decays to 
dilepton final states. In each such decay, one lepton was 
matched to triggering and reconstruction requirements 
so that the other lepton could be used for unbiased effi- 
ciency measurements. This method, known as "tag and 
probe" was used to perform most of the high px lepton 
efficiency measurements used in the dilepton analyses. 
Hadronic jet triggers were studied using events passing 
either muon-based or electron-based triggers. Electron- 
triggered events were required to have exactly one elec- 
tron that was both matched to electron trigger objects 
at all levels and separate from the jet considered in the 
efficiency measurement. 

Efficiency measurements were parameterized in terms 
of the kinematic variables px, f], and <j) of offline recon- 
structed objects. Uncertainties in these parameteriza- 
tions were derived from fits to the observed variable dis- 
tributions. 

Separate efficiencies were estimated for Level I (LI), 
Level 2 (L2), and Level 3 (L3) conditions and the total 
event probability P(L1, L2, L3) was estimated as 

P(L1, L2, L3) = P(L1) ■ P(L2\L1) ■ P(L3\L1, L2), (1) 

where P(L2\L1) and P(L3\L1, L2) are the conditional 
probabilities for an event to satisfy a set of criteria pro- 
vided it has already passed offline selection and previous 
levels of triggering. The overall trigger efficiency for tt 
events was then computed as the luminosity- weighted av- 
erage of the event probabilities associated with each data 
taking period. 



V. EVENT RECONSTRUCTION 
A. Track Reconstruction 

Charged particle trajectories were reconstructed from 
the patterns of energy deposits (or "hits" ) that they left 
in the tracking detectors. Track reconstruction at DO 
involved two distinct steps: track finding and track fit- 
ting. Two complementary track finding algorithms were 
used in event reconstruction. The first is a histogram- 
ming approach based upon the Hough Transform — a 
method originally developed for finding tracks in bub- 
ble chambers [22| . An alternate track-finding approach 
began with groups of hits in the SMT barrels. These 
were fitted to a track hypothesis and the result was used 
to form a road in which to search for hits in additional 
detector layers. 

The candidate track lists resulting from the two track- 
finding approaches were combined and passed to a 
Kalman [23| track fitter. This made use of an interact- 
ing propagator which propagates tracks through the DO 



tracking system while taking into account magnetic cur- 
vature and interactions with detector material. The fitter 
incrementally adds hits to tracks using the input candi- 
dates to define roads. The resulting track fit allows for 
the calculation of optimal track parameters, with errors, 
on any surface. 

The tracking momentum scale was determined by 
comparing the dimuon invariant mass distribution for 
Z — > fifi decays in data with expectation from simulation 
based upon the world average Z boson mass computed by 
the Particle Data Group [3j. In order for the simulated 
track momentum resolutions to match those observed in 
the data sample, an additional random smearing of track 
parameters was performed. 

The measured transverse momentum resolution can be 
expressed as 

<t(1/pt) / (0.003 P t) 2 (0.026) 2 
1/pt V i 4 L-sin0' 1 ' 

Here px is measured in GeV and L is the normalized 
track bending lever arm. L is equal to 1 for tracks with 
\rj\ < 1.62 and is computed as tanfl/tanfl' otherwise (9' 
is the angle at which the track exits the tracker). 



B. Primary Vertex Identification 

The principal task of primary vertex (PV) finding is 
to identify tracks originating from the hard scatter and 
to separate these from tracks generated in superimposed 
minimum bias events. The algorithm first reconstructed 
one or more vertices and then selected the hard scatter 
vertex from among them by considering the px distribu- 
tion and number of tracks associated with each vertex. 

The vertex reconstruction algorithm included three 
steps: track clustering, track selection, and vertex fit- 
ting. First, tracks were clustered in z by considering 
their relative separations. Second, tracks kept for fitting 
were required to have at least 2 SMT hits and px > 0.5 
GeV. Each track's distance of closest approach in the x- 
y plane (dcA) to the nominal interaction position was 
also considered: the d C A significance (S = d C A/o- dcA ) 
for a candidate track had to be less than 3. Finally, for 
every z cluster, an iterative vertex search yielded a ver- 
tex position. A probability that a vertex originated from 
a minimum-bias interaction was assigned based on the 
transverse momenta of its associated tracks. The vertex 
with the lowest probability was selected as the primary, 
or hard scatter, vertex. To further ensure the quality 
of selected primary vertex candidates, we required that 
they be within the SMT fiducial region flzpvi < 60 cm) 
and that they have at least three associated tracks. 

In multijet data events, the position resolution of the 
primary vertex in the transverse plane is around 35 fiin, 
convoluted with a typical beam spot size of 30 [atxi. The 
vertex resolution in the direction of the beam line is ~1 
mm. 
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C. Muon Identification 

Muon identification was based on matches between 
charged particles found in the central tracking system 
and trajectories reconstructed in the muon systems. 
Tracks in the muon detectors comprised straight-line seg- 
ments, or stubs, formed from combinations of scintilla- 
tor and wire chamber hits in a single layer. Stubs were 
formed separately inside and outside the toroid iron and 
then paired together (provided their combinations were 
consistent with expectations for muons originating in the 
interaction region). Pairs were fitted to trajectories us- 
ing knowledge of the toroidal magnetic field and the ex- 
pected effects of energy loss and multiple scattering. The 
resulting local muon momenta were used, along with the 
directions of the muons at the inner surface of the muon 
system, to search for consistent central tracks with which 
to form a global match. Stubs that were not used in form- 
ing local muons were also used in global matches. In all 
cases, the results of the original central track fits were 
taken as the best estimates of muons' momenta, since the 
resolution of the central tracker is far superior to that of 
the local muon system. 

To reduce the impact of muon detector noise, require- 
ments were made on the number and location of scintil- 
lator and wire chamber hits. Two sets of muon quality 
requirements were used in the analyses discussed in this 
paper: tight and loose. Tight muons were required to 
have wire and scintillator hits both inside and outside 
the toroid iron. The loose criteria also accepted muons 
formed from single stubs with both types of hits either 
inside or outside the toroid. 

Central tracks pointing into the fiducial volume of the 
muon detectors (i.e., with |ry| <2) were considered as 
candidates for matches to muon tracks. To ensure that 
a central track was well-reconstructed, the \ 2 P er degree 
of freedom of the Kalman fit used by the central track- 
ing algorithm was required to be less than 4. Consistency 
between candidate tracks and the primary vertex was en- 
sured by two additional cuts: the dcA significance of each 
track must have been less than 3 and the smallest dis- 
tance in z between it and the primary vertex (PV) must 
have been less than 1 cm. The quality criteria applied to 
central tracks and matching local muons are summarized 
in Table [J 

Muons produced in top quark decays can be distin- 
guished from those originating in heavy quark or other 
hadronic decays by a lack of nearby activity in the tracker 
and the calorimeter. This feature motivated two isolation 
criteria used to select signal candidate muons. These in- 
volved summing visible energies over a region around the 
central track associated with the muon. One variable 
was computed by summing the energies of reconstructed 
tracks and the other variable was derived from energy de- 
posited in the calorimeter. For background muons, the 
size of either of these sums is correlated with the muon 
energy, while for signal it is not. Therefore, scaling the 
sums by the pt of the candidate muon generates vari- 



ables that tend to be higher for background than signal. 
A cut on cither of these variables translates to an upper 
limit on surrounding visible energy that increases with 
muon pt- Hence these variables provide more efficient 
criteria than the visible energy sums alone. 
The track-based variable was computed as 

ctrk _ 1 \ ^ trk /o\ 

C halo ~ ~Jt ' 2—1 Pt ' { ' 

Pt A7J<0.5 

Here A1Z was defined for each muon-track pair as 
their separation in r]-<p space, \J Arf + A<fi 2 . The track 
matched to the muon was excluded from the sum. Simi- 
larly, the calorimeter-based isolation was defined as 

££lo = i- E (4) 

P T 0.1<A7J.<0.4 

where the sum was over individual calorimeter cells. In 
each analysis channel that includes muons in its final 
state, signal candidate muons were required to have val- 
ues of both ££* Q and £^f lo that are less than 0.12. This 
requirement was found to reject more than 99% of muons 
originating in hadronic jet decays and to be about 87% 
efficient for a muon coming from a top quark decay. 

D. Electron Identification 

High px electrons were identified by the presence of lo- 
calized energy deposits in the electromagnetic calorime- 
ter. Electron reconstruction began with the formation of 
clusters through the use of a simple cone algorithm that 
grouped calorimeter cells around seed cells having Ex > 
0.5 GeV. Any resulting cluster having Et > 1 GeV was 
then grouped with all EM towers within a cone of ra- 
dius A1Z = 0.4. The centroid of a cluster was calculated 
as the energy weighted mean value of the coordinates of 
its cells in the third layer of the EM calorimeter. Addi- 
tional quality criteria were then applied to reject clusters 
resulting from photons and hadronic activity. 

Electrons (and photons) deposit almost all of their en- 
ergy in the EM section of the calorimeter while hadrons 
typically penetrate into the hadronic sections. Hence an 
electron is expected to have a large EM fraction (Jem), 
which is defined as the ratio of summed energies de- 
posited in the EM layers to the total energy deposited 
inside the clustering cone. 

The longitudinal and lateral shower profiles of an EM 
cluster were required to be compatible with those of 
an electron (or photon). This was done by forming a 
Xcai based on a comparison of the energy depositions in 
each layer of the EM calorimeter and the total energy of 
the shower to average distributions obtained from Monte 
Carlo simulations. 

Electrons produced in top quark decays can be distin- 
guished from those originating in heavy quark or other 
hadronic decays by a lack of nearby activity in the 
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TABLE I: The quality criteria applied to muon candidates. Variable definitions are provided in the text. 



Cut Level 


Requirements 


Loose 


Local muon stubs inside and/or outside toriod iron, 




central track with XKalm»n/U°./. < 4, (y dcA /d CA < 3, 




and Az (track, PV) < 1 cm 


Tight 


Loose and local muon stubs BOTH inside and outside toroid iron 



calorimeter. The electromagnetic isolation fraction (/i so ) 
was used to quantify the degree of isolation of an EM 
cluster and was defined as 

, _ gtot(Ag < 0.4) - E EM {ATZ < 0.2) 
/iso ~ E EM (ATZ<0.2) ' [b) 

where E to t(AlZ < 0.4) is the total energy within a cone 
of size A7Z = 0.4 around the direction of the cluster, and 
Eem(A1Z < 0.2) is the energy in a cone of size A7Z = 0.2 
summed over EM layers only. 

In order to suppress photons and some hadronic jet 
backgrounds, an electron candidate was required to have 
an associated track in the central tracking system within 
\ArjEM,tik\ < 0.05 and | A^M,trk| < 0.05 of the center 
of the EM cluster. 

To further isolate real electrons, an electron likelihood 
formed from seven variables was computed. These vari- 
ables included /em, XZah the ra tio of calorimeter trans- 
verse energy to track transverse momentum (£^ al /p!£ k ), 
the quality of the spatial matching between the central 
track and the EM cluster (XspatialBAf-trk)' tne ^CA of the 
track to the primary vertex, the number of tracks in a 
AR = 0.05 cone, and the sum of the transverse momenta 
of all tracks in a cone of size AR = 0.4 around the EM- 
associated track. Smoothed, normalized distributions of 
each of these variables were made from signal- like (i.e., 
Z —>■ ee) events and background (i.e., QCD dijet) data 
samples. For each discriminating variable Xi, these distri- 
butions provided probabilities P* is {xi) and P^ kg {xi) for 
an EM object to be from a real and a fake electron, 
respectively. The following likelihood discriminant was 
used to distinguish between real electrons and fakes from 
hadronic objects. 

P sig (x) + P b kg(x) 

where x is the vector of likelihood variables. The prob- 
abilities were formed without regard to correlations be- 
tween the likelihood variables, i.e., 

7 

Pi g /bkg(x) = n^s/bkg^)- (j) 

Three classes of electrons were considered on the basis 
of the aforementioned quantities: loose (or EM cluster), 
medium, and tight. The criteria applied to each category 
are listed in Table [TTJ Less than 0.5% of all hadronic 
jets that passed loose EM object identification criteria 



survived the tight cuts. The efficiency of the medium 
quality cuts in data was found to be about 90% in the CC 
and about 63% in the EC. With respect to the medium 
requirements, the additional likelihood cut in the tight 
criteria is about 86% efficient in the CC and 84% efficient 
in the EC. 

The ee analysis selected events with two tight electrons 
and the e + track analysis used one tight electron. The 
e/x analysis had smaller backgrounds and therefore ap- 
plied medium criteria to the signal electron candidate in 
each event. Because of the poor energy resolution for 
electrons reconstructed in the regions between the CC 
and EC sections, all electron-based analyses eliminated 
candidates having 1.1 < \r/\ < 1.5. Electrons with \r/\ > 
2.5 were also removed from consideration in order to sup- 
press multiple scattering backgrounds. 

The EM energy scale was established by requiring that 
the Z boson mass reconstructed in track-matched dielec- 
tron events match the world average Z boson mass com- 
puted by the Particle Data Group [3j . By requiring that 
both electrons in Z candidate events be in either the 
CC or the EC, independent absolute energy scale fac- 
tors were obtained for each portion of the calorimeter. 
These were applied to high px electrons and photons. 
The calibration at lower energies was also checked us- 
ing J/ip —>■ ee decays. The energy resolution for elec- 
trons in the CC or EC is a{E)/E w {15/ VE ® 4)% or 
a(E)/E w (21/s/E © 4)%, respectively (here E is mea- 
sured in units of GeV). 

E. Jet Identification 

Particle jets were reconstructed from energy deposition 
in the calorimeter using a seed-based, improved legacy 
cone algorithm [24j with a cone radius of AR = 0.5. In 
this scheme, seeds were formed by clustering calorimeter 
cells above an energy threshold of 0.5 GeV. All result- 
ing pre-clusters having summed energies above 1.0 GeV 
were then fed into an iterative clustering algorithm. If 
any of the resulting proto-jets shared energy, they were 
either split or merged so that each calorimeter tower was 
assigned to at most one reconstructed jet. Finally, only 
those jets having energies of at least 8 GeV were retained 
for further consideration. 

Additional quality criteria were applied to clustered 
jets to suppress backgrounds originating from noise and 
other instrumental effects. Some of these selection cuts 
were based on variables that discriminate against particu- 
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TABLE II: The quality criteria applied to electron candidates. Variable definitions are provided in the text. 



Cut Level 


Requirements 


Loose (EM cluster) 


^cluster > l£ Qgy^ > Q _ gj ^ < q ^1 < L1 QR L5 < ^1 < 2 5 


Medium 


Loose and /i so < 0.15, Xcai < 50, and central track match 


Tight 


Medium and £ e > 0.85 



lar sources of noise. The coarse hadronic fraction (fen) 
is the fraction of the total energy in a jet that is con- 
tained in the outer, noisier layers of the calorimeter. The 
hot fraction (/hot) is the ratio of the energy of the most 
energetic cell in a jet to that of its next-to-highest energy 
cell. Both /hot and TVgo (the number of cells containing 
90% of the total energy in a jet) were used to suppress 
jets clustered around single cells that fired erroneously. 
Since noise generally did not appear simultaneously in 
the precision readout chain and in the separate Level 1 
trigger readout, the ratio of the Level 1 energy to the 
precision readout energy in a jet (/liv;e t ) is another 
powerful discriminant against jets due to noise. 

Other requirements were made on jet candidates to 
remove clusters that don't originate from partons gen- 
erated in the hard scatter. The EM fraction (/em> see 
Sec. IV Dp . was used to remove reconstructed electrons 
and photons. To eliminate backgrounds from low energy 
multiple interactions, far forward candidates were also 
eliminated. The particular values of all jet quality cuts 
are summarized in Table IIIII 

After these initial cuts, some electrons still remained 
among the reconstructed jet objects. In order to 
avoid the resulting ambiguity, jet candidates overlapping 
medium quality electrons (see Table |LT| within A1Z = 0.5 
were considered only as EM objects. 

A data-to-Monte-Carlo correction factor that accounts 
for possible differences in the jet reconstruction and iden- 
tification efficiencies was determined with back-to-back 
7+jet events by requiring Et balance between the pho- 
ton and the jet. An i/p-dependent scale factor was then 
obtained separately for the CC, EC and ICD regions and 
applied to Monte Carlo. 

A number of effects — including non-linearities in 
calorimeter response, non-instrumented material, and 
noise — can cause the measured jet energy to differ from 
the original particle-level energy. Jet energy scale (JES) 
corrections were applied to adjust jet energies to the par- 
ticle level. Transverse momentum conservation in sam- 
ples of 7 + jet events was used to calibrate JES correc- 
tions in data and simulation. A more detailed description 
of this procedure is available in Ref. [25j . The relative un- 
certainty on the jet energy calibration is rs 7% for jets 
with 20 < p T < 250 GeV. 

Jet momentum resolutions were measured using di- 
jet and 7+jet events. For 50 GeV jets in the CC or 
EC, the resolution was found to be ct(pt)/pt ~ 13% or 
o~(pt)/pt ~ 12%, respectively. 



F. Missing Transverse Energy 

The presence of one or more neutrinos in an event is 
indicated by an imbalance of the visible momentum in 
the transverse plane. Calculation of this quantity be- 
gan with the vector sum of the transverse energies of 
all calorimeter cells surviving various noise suppression 
algorithms, with the possible exception of cells in the 
coarse hadronic layers. These were included only if they 
are clustered within a reconstructed jet. The vector op- 
posite to the result is called the raw missing transverse 
energy ($ T raw )- 

As EM and jet energy scale corrections were applied 
to calorimeter objects, Bt t&w was adjusted through vec- 
tor subtraction. Only jets that pass the quality criteria 
listed in Table [TTTI were used for the hadronic part of this 
correction. The result is called the calorimeter missing 
transverse energy (^?r cal ). 

Since muons are minimum ionizing particles, they de- 
posit only a small fraction of their energy in the calorime- 
ter and Et C31 does not properly account for their pres- 
ence. Therefore, the momenta of all the identified muons 
of tight quality (see Table [T| in an event were subtracted 
vectorially from^r cal after first deducting the muons' ex- 
pected energy depositions in the calorimeter. For the fifj, 
analysis, isolated loose muons which did not pass tight 
quality requirements were also removed from ^r cal . A 
similar procedure was followed for the track identified 
with a signal muon in the /i+track analysis. The fully 
corrected imbalance is simply called the missing trans- 
verse energy ($t)- 

G. b Jet Tagging 

Bottom quark jets were identified using a secondary 
vertex tagging (SVT) algorithm that exploits the long 
lifetime of b hadrons. The algorithm used is the same as 
that used in previouslypublished DO ti production cross 
section measurements d, Q . 

The SVT procedure began by clustering tracks in z 
into track jets. Track jets were reconstructed using a 
A1Z = 0.5 cone algorithm to cluster tracks with px > 0.5 
GeV, at least two SMT hits, and \d C A\ < 0.15 cm. The z- 
projection of a candidate track's dcA onto the beam line 
(zdea) was required to be within 0.4 cm of the z position 
of the PV. 

Within each track jet, tracks having dcA significances 
greater than 3.5, x 2 /d.o.f. less than 3, and transverse 
momentum greater than 1 GeV were paired to form seed 
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TABLE III: The quality criteria applied to jet candidates. Variable definitions are provided in the text. 



Cut 

0.05 < Jem < 0.95 
fen < 0.4 

/hot < 10 

N 90 > 1 

/li E b t >0.4 (H<0.7) 
/liEEt> - 2 (0.7 < M < 1.6) 
M < 2.5 



Target 



Noise and EM particles 
Noise in coarse hadronic layers 
Jets clustered around single cell 
Jets clustered around single cell 
Noise in readout 
Noise in readout 

Extra soft scattering interactions in an event 



vertices. Vertices consistent with having come from 7 
conversions or Kg or A decays were removed from con- 
sideration. Additional tracks pointing to a surviving seed 
were attached to it based on their contribution to the ver- 
tex x 2 - Vertices resulting from this process were selected 
as secondary vertex (SV) candidates based upon the 
collinearity of their component tracks, the \ 2 of their fits, 
and their decay length significances (L xy /<jL my )- Here 
L xy , the decay length, is the distance between the pri- 
mary and secondary vertices in the plane transverse to 
the beam line and ctl xh includes the uncertainty in the 
primary vertex position. The decay length can be posi- 
tive or negative, depending on the sign of its projection 
onto the track jet axis. Secondary vertices correspond- 
ing to the decay of b and, to some extent, c hadrons are 
expected to have large positive decay lengths. 

A calorimeter jet was tagged as a & jet if a secondary 
vertex with L xy /o~L xy > 7 was found within A1Z < 0.5. 
If a jet contained at least one secondary vertex with 
L xy /o~L xy < —7, the jet was labeled negatively tagged. 
Negative tags resulted from fake or mis-reconstructed 
tracks, or from the effects of multiple scattering in de- 
tector material. Negative tags were used to estimate the 
probability to misidentify a light flavor (a u, d, or s quark 
or a gluon) jet as a b jet (the mis-tagging rate). 

The overall event tagging probability for a particu- 
lar process depends upon the flavor composition of the 
jets in the final state and on the event kinematics. This 
probability was estimated through the application of tag- 
ging rates measured in data to each jet in simulation. A 
brief description of tagging probability measurements is 
given here, and a more detailed discussion is available in 
Ref. @. 

In order to decouple tagging efficiency measurements 
from detector geometry effects, jets considered for tag- 
ging were required to pass additional taggability criteria. 
A taggable jet had to have its axis matched to within 
ATI < 0.5 with the axis of a track jet. The SMT hit 
requirement for tracks in track jets means that most tag- 
gable jets are associated with PV z positions within «36 
cm of the center of DO detector. Within this range, jets 
with momenta above 30 GeV typically have taggabilities 
of greater than 90%. 

The b tagging efficiency was estimated using a sample 
of dijet events enriched in semileptonic decays of bottom 
and charm hadrons by the requirement that one jet have 



an associated muon. The heavy flavor content of this 
sample was further enriched by requiring that the op- 
posite jet was tagged, either with the SVT algorithm or 
with a soft lepton tag. Soft lepton tagging requires that a 
muon with a momentum component along the direction 
orthogonal to the jet axis of at least 0.7 GeV be found 
within A1Z = 0.5. In order to extract the SVT tagging ef- 
ficiency, both tagging algorithms were applied separately 
and together to the dijet sample. The resulting b tagging 
efficiency depends upon both rj and px, and a typical 40 
GeV taggable jet from top quark decay is tagged about 
40% of the time. Charm quark jets have tagging efficien- 
cies around 20% as large as those for b quarks. Light 
quark mis-tagging probabilities are on the order of 0.1%. 

VI. ANALYSES 
A. General Considerations 

As discussed in Sec. II B[ ti decays to dilepton final 
states are characterized by the presence of two high pt 
charged leptons, significant $t from two neutrinos, and 
two or more jets (from b quark fragmentation and ini- 
tial and final state radiation). Analysis of all four chan- 
nels described here therefore begins by requiring that sig- 
natures consistent with two isolated, oppositely charged 
leptons and at least two jets be reconstructed. 

The backgrounds at this stage vary with the channel 
considered, but generally include Drell-Yan production 
of (Z/7*)+jets, diboson production (i.e., WW or WZ) 
with jets, and leptonic VF+jets events in which another 
lepton arises from the misidentification of one of the jets. 
Resonant production of Z bosons which decay into elec- 
tron or muon pairs is the dominant background for the 
ee and fifi channels, each of which employed cuts based 
upon $t and the invariant mass of the lepton pair to 
target this process. 

The remaining backgrounds were removed using a com- 
bination of kinematic and topological constraints. These 
include a summed transverse energy Ht, defined as 

2 

H T =t$+Y,T%\ (8) 

»=i 

where l\ denotes the highest pt lepton, and i in the sum- 
mation extends over the two highest pt jets in the event. 
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For the e/i analysis, a cut on Ht was found to be more 
effective than one onj)?T hi rejecting Z — ► rr background. 

Analyses using fully reconstructed leptons (i.e., the 
e/i, ee, and /i/i channels) were optimized separately to 
achieve the best possible performance using kinematic 
and topological quantities to supress dominant back- 
grounds. In order to recover some of the efficiency lost 
due to lepton identification requirements, an alternate 
approach using one fully reconstructed lepton and one 
central track was taken. To contend with the additional 
background let in by the lack of lepton identification re- 
quirements on the second lepton, the ^+track analysis 
required that at least one jet pass explicit b quark tag- 
ging requirements. The ee, fifi, and i?-|-track selections 
are not completely orthogonal, and their overlaps were 
accounted for in the combined cross section calculation 
discussed in Sec. I VIII The low-background e/i analysis 
was designed to have no overlap with the other channels. 



B. The e/i Channel 

The signature for an e/i event consists of one high pt 
isolated electron, one high px isolated muon, and two or 
more jets. The major backgrounds to this channel come 
from Drell-Yan production of r pairs which in turn de- 
cay to produce an e/i pair (i.e., Z/j* — > rr — > e/i + X) 
and WW production with jets. There are additional 
backgrounds present from misidentificd leptons, particu- 
larly electrons. These are mostly W(— ■> /zz/)+3 jet events 
in which one of the jets was misidentified an electron. 
Hereafter, objects misidentified as electrons will also be 
referred to as fake electrons. 

Offline selection began with medium electron and tight 
muon identification cuts (see Sees. IV Dl and IV CI for de- 
tails). To reject bremsstrahlung events, in which the 
muon emits a photon that is mistakenly identified as an 
electron, the candidate electron and muon were not al- 
lowed to share a common track in the central tracking 
detectors. Additionally, the candidate electron and muon 
were required to be matched to tracks of opposite charge. 

After this initial lepton selection, the sample was domi- 
nated by background consisting of roughly equal amounts 
of misidentified leptons and physics processes leading to 
legitimate e/i pairs. The backgrounds generally contain 
jets arising from initial and final state radiation. These 
tend to be softer in pt than the b jets that are generated 
in it decays. Requiring that two or more jets (Sec. IV E[) 
be found with pr of at least 20 GeV reduces both back- 
grounds by more than a factor of 50 while preserving 
more than two-thirds of the signal. 

In addition to basic lepton and jet identification and 
energy requirements, the use of event-wide selection crite- 
ria was found to improve the expected significance of the 
result. Several variables were considered, including 
Ht, and the transverse mass of the combination of the 
leptons and the $t- The performances of cuts on these 
quantities, whether alone or in combination, were eval- 
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FIG. 1: Ht distributions in the e/i channel for expected sig- 
nal (assuming mt = 175 GeV and er t j = 7 pb), expected 
background, and data after both lepton and jet requirements 
(corresponding to the second line of Table HV)) . 



uated using the expected significance of the background 
subtracted yield, including both the statistical uncertain- 
ties of the yields and a term reflecting the dominant jet 
energy scale systematic uncertainty in the total back- 
ground estimate. It was found that requiring at least 122 
GeV of Ht gave the best performance. Figure Q] shows 
the Ht distributions for signal and background before 
this cut and illustrates its ability to discriminate signal 
from background. Table []V] shows the impact of the Ht 
cut on expected signal and background yields. 

After all cuts are applied, 21 events remain in the data. 
Table IVl shows the expected background and signal (as- 
suming TO f = 175 GeV and cr t j = 7 pb) contributions to 
the final sample. 

Contributions from the processes Z — > tt — ► e/i + X 
and WW/WZ — >■ efi+X were calculated using the Monte 
Carlo samples discussed in Sec. Mil The background from 
fake electrons was estimated by performing an extended 
unbinned likelihood fit to the observed electron likelihood 
(Sec. rVD|) distributio n in events passing all selection cri- 
teria. The distribution, which is shown in Fig. [21 was 
fitted using a likelihood given by 



JV 



£ = JJ^e-SC^i) + n i3ke B(xi)]- 



i=i 



-(n c +n fako ) 



(9) 



where i is an index that runs over all selected events, 
Xi is the corresponding observed value of the electron 
likelihood, N is the total number of events, n e is the 
number of events with signal- like electrons, nf a k e is the 
number of events having fake electrons, and S and B 
are the signal and background probability distribution 
functions, respectively. The event counts n e and nf a k e 
were allowed to float in the fit. 
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TABLE IV: Numbers of observed and expected e/j, events passing the analysis cuts. The instrumental background is from fake 
electrons. Expected number of tt events are for m t = 175 GeV and a t f = 7 pb. Uncertainties correspond to statistical and 
systematic contributions added in quadrature. 







Total 


Instrumental 


Physics 






Data 


sig+bkg 


bkg 


bkg 


tt 


Trigger, 1 e, p T > 15 GeV 












+ >1 At, Pt > 15 GeV 
+ >2 jets, p& > 20 GeV 
+ H T > 122 GeV 


24 
21 


22.0 + 2g 
17.7t*5 


o n+2.8 
°- z -2.0 
1+2.5 
z - 1 -1.7 


2 5 +0 ' 7 


13.8+i'? 

i3.iii-;s 



TABLE V: A more detailed listing of the expected e/j, signal 
and background yields presented on the last line of Table IIVI 
The expected number of tt events is calculated assuming mt = 
175 GeV and a t t = 7 pb. Uncertainties include statistical and 
systematic contributions added in quadrature. 



Process 



Expected number of e/i + X events 



a (mc) 

Z -> tt -> e/i + X (MC) 
WW/WZ -^en + X (MC) 
Fake leptons (data) 
Total background 



13.09™ 

+0.40 
0.42 
2.50 
1.66 
+2.56 
1.77 



0.99 
2.14 
4.58 



probability distribution function was determined using 
events passing all signal selection criteria but the jet re- 
quirements and using an anti-isolation cut on the muon. 
The contribution from signal-like electrons to the re- 
sulting sample was found to be less then 0.5%. Note 
that the fake electron estimate resulting from Eq. [5] in- 
cludes backgrounds containing both real and fake isolated 
muons. The contribution from events containing legiti- 
mately identified electrons and falsely isolated muons was 
also investigated and found to be negligible. 



The ee Channel 
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FIG. 2: The electron likelihood distribution in the e/x channel. 
Data passing all selection criteria are shown with points. The 
fake electron background estimated by a fit to the data points 
is shown, along with the simulated likelihood distributions 
for the each of the estimated signal and other background 
contributions, as shaded histograms. 



The probability distribution functions used in the fit 
were determined with separate data samples enhanced 
in signal-like and background-like electrons. The signal 
probability distribution function came from a fit to the 
likelihood distribution of electrons in oppositely signed 
dielectron events selected using standard electron iden- 
tification cuts and having low $t- The background 



The signature for an ee event consists of two high pr 
isolated electrons, at least two high px jets, and substan- 
tial The main background to this signature arises 
from Drell-Yan production of dielectrons (Z/j* — > ee). 
Although this process produces no real ~$t-, mismeasured 
l^T can originate from misreconstructed jet or electron 
energies or from noise in the calorimeter. Other back- 
grounds include Drell-Yan production of r pairs which 
further decay to dielectrons (Z/j* — > tt — > ee + X) 
as well as diboson (WW/WZ) production associated 
with jets. There are also small backgrounds from W(—> 
e^)+multijet and QCD multijet events in which one or 
two jets are misidentified as isolated electrons. The back- 
ground from heavy flavor (cc, bb) production is negligible 
since electrons from these decays are typically soft and 
non-isolated. 

Event selection began with two tight electrons, as de- 
scribed in Sec. IV D[ each having px > 15 GeV. The elec- 
trons were also required to have matching central tracks 
of opposite charge. This initial selection essentially elim- 
inated any background from heavy flavor production and 
significantly reduced the background from misidentified 
electrons. The additional requirement that two jets be 
found, each with px > 20 GeV, generated a sample domi- 
nated by Drell-Yan Z/j* production with associated jets. 
Table IVII shows the ee sample composition at this and 
subsequent stages of selection. 

Simultaneous cuts on the $t and the dielectron invari- 
ant mass (M ee ) provide a powerful way to suppress most 
of the Z/"f* — > ee background. We vetoed events with 
M ee values near the Z boson mass (i.e., 80 GeV < M ee < 
100 GeV) and required # T > 35 GeV (fi T > 40 GeV) for 
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TABLE VI: Numbers of observed and expected ee events passing the analysis cuts. The instrumental background includes 
events containing misidentified electrons and misreconstructed $t- Expected number of tt events are for m t = 175 GeV and 
a t i = 7 pb. Uncertainties correspond to statistical and systematic contributions added in quadrature. 





Data 


Total 
sig+bkg 


Instrumental 
bkg 


Physics 
bkg 


tt 


Trigger, N e >2,p e T > 15 GeV 

+ >2 jets, p£* > 20 GeV 
+ M ee cut 
+ $ T cut 
+ sphericity cut 


369 
88 
5 
5 


428.4+™-? 
106.2+^ 

r y+0.5 
c 9+0-5 


415.9+™-? 
98.61^:° 
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0.5 ±0.2 
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i-9±g:g 
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c 7+0.5 
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FIG. 3: Expected and observed |?t distributions in dielectron events before (left) and after (right) cutting on M ee and $t- 



M ee > 100 GeV (M ee < 80 GeV). These requirements 
effectively suppressed the Z/"f* — > ee background and 
brought other backgrounds to a manageable level while 
preserving significant signal efficiency. Their effect is il- 
lustrated in Fig. [3] 

Finally, a cut on the sphericity of the event was applied 
in order to take advantage of the topological peculiarities 
of tt events and gain even more discrimination between 
signal and background. Sphericity (S) is defined as 



(10) 



where e\ and £2 are the two leading eigenvalues of the 
event- normalized momentum tensor [26j. The tensor 
(■M-xy) is calculated as 



M xy = 



(ii) 



where the index i runs over the leading two electrons and 
the leading two jets in the event. 

Sphericity can take values between and 1. The ap- 
plied cut of S > 0.15 rejects events in which jets are 
produced in a planar geometry due to gluon radiation 



TABLE VII: A more detailed listing of the expected ee signal 
and background yields presented on the last line of Table IVII 
The expected number of tt events is calculated assuming mt = 
175 GeV and a t i = 7 pb. Uncertainties include statistical and 
systematic contributions added in quadrature. 



Process 



Expected number of ee + X events 



tt (MC) 

Z -> tt -> ee + X (MC) 
WW/WZ —* ee + X (MC) 
Mismeasured fir (data) 
Fake electron (data) 
Total background 



4.04Z 
0.35j 



4+0.40 
-0.46 
S+0.11 
-0.15 

0.45 ±0.15 
0.09 ±0.03 



and provides a reasonable reduction in most of the back- 
grounds. The cut value was chosen using a figure of 
merit related to the expected statistical significance of 
the background subtracted signal. The observed and ex- 
pected sphericity distributions for events passing the M ee 
and l^T cuts are shown in Fig. |4j 

After all cuts, five events remain in the data. Table [VTll 
shows the the expected background and signal (assuming 
m t = 175 GeV and cr t j = 7 pb) contributions to the final 
sample. 
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FIG. 4: The observed and predicted sphericity distributions 
of dielectron events after the M ee and $t cuts. 
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FIG. 5: The probability to pass the J?t selection for events 
with two or more jets in Z/^y* — > e + e~ data, in 7 + 2 jets 
data, and for Z/j* —> e + e~ Monte Carlo. 



The background contributions from Z/j* — > tt and 
WW/W Z{-^ ee + X) decays were estimated using the 
Monte Carlo samples described in Sec. IIIIl Simu- 
lated data were also used to estimate contributions from 
Z/7* — ► ee decays for cut levels prior to the 1£t require- 
ment in Table IVTl The background due to mismeasured 
$t was estimated from data using the observed misrecon- 
struction probability (/met) in 7+2 jets events selected 
to have kinematics and resolutions similar to the Z/j* 
backgrounds. Figure [5] shows that the $t spectrum in 
Z/7* — > ee events with two or more jets is in good 
agreement with that observed in the 7+2 jets sample. 
We therefore estimated /met as the fraction of 7+2 jets 
events passing the $t selection. 



The $t misreconstruction probability was multiplied 
by the number of events that failed the J£t cut but passed 
all other selections. Thus the number of expected $t 
background events summed over the low and high M ee 
regions is 



'"MET — JV M cc <80GeV 



xA 4 



40 GoV , «r „ f35 GoV 

MET +-< v M„>100GeV x /MET • 

(12) 

The background due to electron misidentification was 
also obtained from data by calculating the fake electron 
probability, f e . This was derived from a control sample 
containing two loose EM objects and passing signal di- 
electron triggers. Additional cuts on $t and M ee were 
used to suppress contributions from signal-like electrons, 
and the resulting sample was found to be completely 
dominated by fake electrons. 

The predicted number of fake electrons in the final 
sample (iV e fake ) was obtained by multiplying the number 
of events with one loose and one tight electron by f e . At 
this stage, JVg ako contains both W+jet and QCD multi- 
jet backgrounds. The latter enters into the sample when 
two jets mimic the signal electron signature and misre- 
constructed $t allows the events to pass full selection 
criteria. Since this background was counted along with 
the misreconstructed $t background obtained from the 
data, it was removed from the A^ a ke estimate to avoid 
double counting. This was achieved by loosening iden- 
tification cuts on both electron candidates in the final 
sample and scaling the yield by the square of f e . This 
led to an estimate of Nqcd that was used to correct 

jV fake_ 



D. The /i/x Channel 

The signature for a fj.fi event consists of two high 
Pt isolated muons, two high px jets, and large $t- 
The dominant background in the dimuon channel comes 
from Drell-Yan production of muon pairs (Z/j* — * i-j.fi) 
in which misreconstructed objects give rise to mismea- 
sured $t- Other backgrounds include Drell-Yan pro- 
duction of t pairs which decay to produce muon pairs 
(Z/7* — > tt — > fifi + X) as well as WW and WZ produc- 
tion with jets. There are also backgrounds from multijct 
and VF+jets events, in which one or more muons from 
heavy flavor decay pass isolation criteria. Hereafter, this 
class of muons will also be referred to as fake isolated 
muons. 

Candidate events were required to contain two loose, 
isolated muons (as described in Sec. IV C[) with match- 
ing central tracks of opposite charge and two jets with 
p T > 20 GeV. The first line of Table IVlTTl shows signal 
and background yields after this initial event selection. 
The data selected are heavily dominated by misidcntificd 
Drell-Yan events. Two additional selection requirements 
were designed to specifically target this background. The 
first was a contour cut made in the plane formed by 
event $t and the opening angle in <f> between the leading 
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TABLE VIII: Numbers of observed and expected fifj, events passing the analysis cuts. The instrumental background includes 
events containing fake isolated muons and misreconstructed Expected number of tt events are for m t = 175 GeV and 
a t t — 7 pb. Uncertainties correspond to statistical and systematic contributions added in quadrature. 
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Data 
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+ X 2 >4 
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FIG. 6: Expected and observed $t distributions in dimuon events before (left) and after (right) the two-dimensional cut in the 
A0(/ii ca ding,^T) versus $ T plane. 



Pt muon and the $t- Correlations between these two 
variables are caused by the misreconstruction of central 
tracks matched to muons. An event having $t less than 
45 GeV was immediately rejected. This cut was further 
tightened at low and high values of A0(^i ea ding,^T) to 
$ T > 90 GeV and $ T > 95 GeV respectively. Events 
with A0(^ieading,^T) > 175° were removed. As can be 
seen in Fig. [F5] and the second line of Table IVIII1 the 
contour cut effectively suppressed the misidentified back- 
ground. 

Further background rejection was achieved by cutting 
on the compatibility of an event with the Z — * hy- 
pothesis. To this end, a \ 2 was formed using a Z boson 
mass constraint and the measured muon momentum res- 
olution. The resulting variable, shown in Fig. [71 accounts 
for the pt and r\ dependence of the tracking resolution 
and was found to perform better than a simple dimuon 
mass cut. The final cut value (\ 2 > 4) and the location 
and shape of the contour cut described above were cho- 
sen using a grid search over Monte Carlo predictions of 
signal and background yields. Both cuts were varied si- 
multaneously in the search and the best combination was 
chosen using a figure of merit related to the expected sta- 
tistical significance of the background subtracted signal 
and including the expected uncertainty on the Z + 2 jets 
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FIG. 7: The observed and predicted \ 2 distributions of 
dimuon events before the \ 2 cu t- 



background prediction. 

After all cuts, two events remain in the data. Ta- 
ble IIXI shows the expected background and signal (as- 
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TABLE IX: A more detailed listing of the expected fifj, signal 
and background yields presented on the last line of Table IVTTT1 
The expected number of ti events was calculated assuming 
m t = 175 GeV and a t i = 7 pb. Uncertainties include statis- 
tical and systematic contributions added in quadrature. 



Process 



Expected number of /x/i + X events 



ti (MC) 

WW/WZ -> nfi + X (MC) 
Z -> hh/tt(-> w)+X (MC) 
Isolation fakes (data) 
Total bkg 



2.96 4 



FDT3T 
-0.35 

u - iy -0.07 

u -*'-0.18 
,+0.01 
0.01 

+0.24 



0.01 
0.67 



suming m t = 175 GeV and cr t j = 7 pb) contributions to 
the final sample. 

Background estimates for Z — ► /i/i, Z — > tt — ► + A", 
and WW/WZ production were calculated using the 
Monte Carlo samples discussed in Sec. IIII1 The back- 
ground from fake isolated muons was estimated using a 
procedure requiring two samples of events: a "tight" sam- 
ple containing Nt events passing all the dimuon selec- 
tion criteria and a "loose" sample of Nl events for which 
only one muon was required to pass isolation cuts. These 
event counts are related to the numbers of events with 
signal- like muons (N s i) and events with background- like 
muons (Nbi) via the relations: 



and 



N T = e sig N st + f^N u . 



(13) 



(14) 



Here e s i g is the probability for signal-like muons to pass 
isolation requirements and is the probability for muons 
in background events to pass isolation requirements. 
Equations [TS] and [TH can be solved for the falsely iso- 
lated muon background in the fully selected sample (i.e., 
UN U ). 

The faking probability was estimated as the isolation 
efficiency for the second highest pt muon in dimuon 
events. To eliminate bias from W and Z boson decays, 
events in the sample were also required to have non- 
isolated leading muons and dimuon masses less than 70 
GeV. The isolation efficiency for signal-like muons was 
estimated with ALPGEN ti — * /i/i + X Monte Carlo. 



E. The ^+track Channel 

The final state in the £+track channel is character- 
ized by two oppositely charged high px leptons, one ex- 
plicitly reconstructed as an electron or a muon and the 
other identified as an isolated track. Requiring an iso- 
lated track as opposed to a fully reconstructed lepton 
allows the recovery of some events lost due to lepton re- 
construction inefficiency and also adds a small number 
of events with a hadronically decaying tau lepton to the 



sample (5%/4% of the e+track//i+track data). Events 
were also required to contain two high pr jets and a large 
$t ■ The dominant background in this channel originates 
from Drcll-Yan production of lepton pairs {Z/"f* — > ee or 
Z/j* — ► fifj,) where misreconstructed objects, resolution 
effects or noise can give rise to mismeasured $t- Addi- 
tional instrumental backgrounds arise from multijet and 
W+jets events, where a jet can be misidentified as an 
isolated lepton or track. Sources of irreducible (physics) 
backgrounds are Z/j* — ► tt and WW/WZ production 
in association with jets. 

Selected events were required to have one tight, iso- 
lated lepton (electron or muon) with px > 15 GeV and 
an oppositely charged isolated track with pt > 15 GeV. 
Neither the lepton nor the isolated track was allowed to 
be found within the cone of a reconstructed jet. In addi- 
tion, the lepton and the isolated track were required to be 
separated in A1Z by requiring A7£(lepton, track) > 0.5. 

A track was considered isolated if £^ a , defined in 
Eq. [3l was less than 0.12. The quality criteria applied 
to the isolated track were identical to those applied to 
the central muon tracks, with the exception of the dcA 
significance requirement, where <Jd CA /dcA < 5 was used. 

The candidate events were further required to have 
two or more jets, each with px > 20 GeV. The sam- 
ple selected by the above requirements is dominated by 
Z/-f* background events. This background was partly 
removed by requirements on the $t hi each event. For 
the /j,+track channel, $t > 35 GeV (> 25 GeV) was 
required for 70 GeV< M^ jtrk < 110 GeV (M M , trk < 70 
GeV or M M , t rk > 110 GeV). The corresponding cuts in 
the e+track channel were $t > 20 GeV (> 15 GeV) 
for 70 GeV < M ejtrk < 100 GeV (M e , trk < 70 GeV 
or M ejtr k > 100 GeV). The tighter requirements for the 
/z+track channel reflect the impact of the muon-matched 
central track px resolution on the reconstructed trans- 
verse energy. 

While the preceding requirements are effective, the 
most powerful discriminant used to suppress backgrounds 
in this analysis is the requirement of at least one b tagged 
jet, since jets in background events originate predomi- 
nantly from light (u, d, or s) quarks or gluons. The b 
tagging algorithm and its performance are discussed in 
Sec. IV Gl Figure [5] shows the numbers of observed and 
predicted events with two or more jets as a function of 
Ht (defined in Eq. [5]) before and after the secondary ver- 
tex requirement. The impact of the tagging requirement 
on the backgrounds is clearly visible. Tables |X] and IXll 
show the impact of the J^t and b tagging requirements 
on the expected signal and background yields. 

After application of all selection criteria, 17 events re- 
main in the data. Tables IXlTl and IXIIII show the expected 
signal and background contributions to the final sample. 

The b quark and c quark tagging efficiencies and mistag 
rate are parameterized as functions of jet px and \rj\ (see 
Sec. IV G|) . The tagging probabilities for ti events were 
estimated by applying these parameterizations to jets in 
simulated events. 
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FIG. 8: J/t spectra of observed and predicted events in the ^+track channel before (left) and after (right) the secondary vertex 
requirement. 



TABLE X: Numbers of observed and expected e+track events passing the analysis cuts. The instrumental background includes 
events containing misidentified electrons and misreconstructed ]£t- Expected numbers of tt events are for mt = 175 GeV and 
au = 7 pb. Uncertainties correspond to statistical and systematic contributions added in quadrature. 





Data 


Total 
sig+bkg 


Instrumental 
bkg 


Physics 
bkg 


tt 


Trigger, e+track, p?/ trk > 15 GeV 
+ >2 jets, p£* > 20 GeV 

+ |?T CUt 

+ > 1 b tagged jet 


436 
85 
11 


442.21™;? 

no n+20.8 
^■^-13.8 

io.9±l:S 


422.6 + ^| 
74.5tf 3 i 

9 7 +0.9 
z - ' -0.7 


4 7 +0.9 
q O+0.7 

1+ 00 
u - 1 -o.o 


15.0± ;° 
I3.8i ;° 

8.it° :S 



The Z/j* — ► ee/iijjL, Z/^j* — + rr, and diboson back- 
grounds were estimated using the simulated samples dis- 
cussed in Sec. IIIII The b tagging probability for Z/j* — > 
ee/fifj, events was estimated using a control sample se- 
lected with all /?+track event selection criteria with the 
exception of the $t cut, which was reversed. As a cross 
check, separate b tagging probabilities were measured for 
the e+track and /^+track channels and were found to be 
consistent. The b tagging efficiency for Z/j* — > tt was 
taken to be the same as for Z/j* — > ee,/z^ events. For 
diboson events, the b tagging efficiency was assumed to 
be the same as for W boson events with associated jet 
production. 

The contribution from W boson and multijet events 
containing fake isolated leptons and/or tracks was esti- 
mated using specially selected data samples. The method 
employed is similar to that described for similar back- 
grounds in the dimuon analysis (Sec. lVIDj) . Because the 
effects of loosening isolation requirements on leptons and 
tracks were included separately, the system of equations 
1131 and [TJ] expands to include four equations. The un- 
knowns in each expression include the number of events 
with fake isolated tracks and fake isolated leptons (N^), 
the event count with fake isolated tracks and real isolated 
leptons (Nyf), the number of events with real isolated 



tracks and fake isolated leptons (iVj?'), and the event 
count with real isolated tracks and real isolated leptons 
(N^f). Once these quantities are obtained, backgrounds 
due to fake isolated leptons and/or tracks can be com- 
puted as 

AWts - Ay? + N$ (15) 

and 

N QC u=N f f l (16) 

In order to limit the impact of statistical fluctuations, 
this procedure was applied to events selected without b 
tagging. The final background estimates were then de- 
rived using separate tagging efficiencies for VF+jets and 
multijet backgrounds. 

Knowledge of the efficiencies of the tight track and lep- 
ton criteria relative to their loose counterparts for both 
signal-like and background-like objects is required to ex- 
tract the unknown yields from observed event counts. 
The efficiency for a signal-like loose electron, muon or 
track to pass each corresponding tight criterion was de- 
termined from simulated samples of Z/j* — > ee and fifi 
events. The efficiencies for loose fake leptons were mea- 
sured in a multijet data sample obtained by selecting 
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TABLE XI: Numbers of observed and expected /i+track events passing the analysis cuts. The instrumental background includes 
events containing fake isolated muons and misreconstructed |?t- Expected numbers of ti events are for m t = 175 GeV and 
o t f = 7 pb. Uncertainties correspond to statistical and systematic contributions added in quadrature. 





Data 


Total 
sig+bkg 


Instrumental 
bkg 


Physics 
bkg 




Trigger, ^i+track, p^ t,k > 15 GeV 

+ >2 jets, p£* > 20 GeV 
+ $ T cut 
+ > 1 b tagged jet 


480 
56 
6 


483.5tf 2 j 
63.8t£i 

o q+0.8 

8-3-0.8 


465.8t 8 8 ^ 

-0. i 
i-9±S:S 


4.5+?-? 
2 6+° :7 

u^-o.o 


I3.it!; , 
10.8±g;g 



TABLE XII: A more detailed listing of the expected e+track 
signal and background yields presented on the last line of 
Table [X] The expected number of ti events is calculated as- 
suming mt = 175 GeV and a t t = 7 pb. Uncertainties include 
statistical and systematic contributions added in quadrature. 



Process 



Expected number of e+track events 



ti (MC) 

WW (MC) 

Z -> ee,fj,fi (MC) 

Z^TT (MC) 

VK/multijet (data) 
Total bkg 



S +U.U8 

0.02+ 109 



2.29 
0.12 
0.42 
2.85 



1.09 
+0.40 
0.31 
+0.31 
0.35 
+0.37 
0.37 
+0.33 
0.27 



TABLE XIII: A more detailed listing of the expected /z+track 
signal and background yields presented on the last line of Ta- 
ble IXII The expected number of ti events is calculated as- 
suming mt = 175 GeV and a t t = 7 pb. Uncertainties include 
statistical and systematic contributions added in quadrature. 



Process 



Expected number of /z+track events 



ti (MC) 

WW (MC) 

Z -> ee, (MC) 

Z^TT (MC) 

W /multijet (data) 
Total bkg 



fi 2 q+U.U9 
D - zy -0.09 



0.01 
1.83 
0.08 
0.08 
2.00 



1.10 
1.12 
+0.31 
0.30 
+0.36 
0.46 
+ 0.88 
0.88 
0.2') 
0.30 



TABLE XIV: Numbers of observed events, expected back- 
ground yields, the product of ti selection efficiency times 
branching ratio, and the integrated luminosity for each anal- 
ysis channel. The branching fractions for the e/z, ee and \i\x 
channels considers the decays ti — > bbWW — > e/i/ee//i/i + X 
respectively. Both e+track and /i+track channels consider 
ti —f bbWW — > It + X decays (£ — e, fi, r) with the r leptons 
decaying both leptonically and hadronically. 



Channel 


N obs 


AW 


e x B (%) 


JCdt (pb" 1 ) 


e/i 


21 




0.44 ± 0.04 


427 ± 26 


ee 


5 


1 1 n+0.22 
1,J -^— 0.27 


0.13 + 0.02 


446 ± 27 


W 


2 


n R7+ - 24 

U.D(_ 2 2 


0.10 + 0.02 


421 ± 26 


e+track 


11 


9 oc+0.33 
Z.OO_ 27 


0.27 + 0.02 


425 ± 26 


/z+track 


6 


2 00+ ' 29 

z - uu -0.30 


0.21 + 0.02 


422 ± 26 



F. Summary 



Table IXTVl presents a summary of event counts ob- 
served in data, expected background yields, the products 
of ti selection efficiencies and branching ratios, and lu- 
minosities for each of the five dilepton analysis channels. 
These quantities enter into the top quark pair produc- 
tion cross section calculations discussed in Sees. I VII Al 
and lVlTCl 



VII. CROSS SECTION CALCULATION 



lepton events in a low $t (<10 GeV) region. Biases 
from Z boson decays were removed by eliminating events 
with two like-flavor leptons or with an additional isolated 
track that, when paired with the lepton, formed an in- 
variant mass consistent with Mz- Efficiencies for loose 
fake tracks were measured in similarly chosen samples. 

Secondary vertex tagging efficiencies for W^+jets events 
were measured using single-lepton events selected with 
the same Z boson rejection criteria used for the mul- 
tijet sample described above. Additional biases from 
the presence of top quark pairs were accounted for by 
subtracting predicted ti contributions calculated with an 
assumed production cross section of 7 pb. The event 
tagging probabilities for multijet events were determined 
using the same multijet data samples used to estimate 
isolation faking probabilities. 



A. Individual Channel Cross Sections 

To estimate the production cross section a 3 for an indi- 
vidual dilepton channel j, the following likelihood func- 
tion was defined: 

{A r ° bs , N^ kfS , Ej x Bj,Cj}) = V(Nf s ,n 3 ) 

N° bB 

TI 

!\7obs| 'V > 

3 ' 

where P(N° hs ,rij) is the Poisson probability to observe 
N° hs events given an expected combined signal and back- 
ground yield of 

nj = a 3 { £j xB 3 )C +N^. (18) 
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TABLE XV: The tt production cross sections at y^i = 
1.96 TeV and for a top quark mass of 175 GeV as measured 
in each analysis channel. The first uncertainty listed for each 
result is statistical in origin. The second uncertainty is the 
combined effect of all systematics, excluding the uncertainty 
on the luminosity measurement. The final error listed is from 
the luminosity measurement. The origins of the systematic 
uncertainties are discussed in Sec. IVH Bl 

Channel au (pb) 

~e/l 8.8 +23 -i't ±0-5 

ee 6.7 tH ±li ±0.4 

w 3.i tii ±8:1 ±0.2 

e+track 7.1 t^e -li ±°- 4 

^i+track 4.5 ±f;J ±° 9 9 ±0.3 



Here Cj is the luminosity, Cj x Bj is the product of se- 
lection efficiency and branching fraction, and A^ bkg is the 
expected background for channel j (see Table IXTVT) . The 
cross section is then extracted by minimizing the negative 
log-likelihood function, 

-logL(a j ,{Nf s ,N^,e j x B h C j }). (19) 
The results are presented in Table IXVI 



B. Systematic Uncertainties 

Systematic uncertainties for the analyses can be 
broadly grouped into those related to signal acceptance 
calculations and those concerning overall background es- 
timates. Brief descriptions of the sources of systematics 
are provided below. 

• Primary vertex identification 

A correction to the simulated efficiency for pri- 
mary vertex selection was estimated by compar- 
ing its value in Z — > ee//x/z decays in data and 
Monte Carlo. In order to quote a systematic uncer- 
tainty related to this correction, the ratio of data 
and Monte Carlo efficiencies was varied by ±ler, 
and signal efficiencies and expected background 
yields were re-computed. The ultimate origin of 
the uncertainty a is the statistical limitations of 
the Z — > ee/ mi data samples. 

• Lepton identification 

For electrons and muons, uncertainties related to 
the identification and selection criteria described 
in Sees. IV CI and IV Dl were estimated using control 
samples of Z — ► ee//i/i decays. The tag and probe 
technique discussed in Sec. |W] was used to com- 
pute the effects of each criterion in data and Monte 
Carlo, and the ratio of the resulting efficiencies was 
used to correct the simulation. When a correction 
was found to depend on object kinematics, it was 



binned appropriately. The corresponding system- 
atic uncertainty was computed by varying the cor- 
rection by ±1(7, where a arose from the statistical 
limitations of the Z boson control samples in data. 

• Track reconstruction 

Analyses using muons and the £+track channels 
include uncertainties associated with central track 
reconstruction. Chief among these is the uncer- 
tainty in the track smearing procedure discussed in 
Sec. IV Al Signal efficiency and background yield 
calculations were repeated with the smearing pa- 
rameters varied according to their uncertainties, 
whose ultimate origin is in the parameterization 
of the smearing functions and the size of the data 
samples used to tune them. Because of the signif- 
icance of Bremsstrahlung energy loss for electrons, 
separate uncertainties were used for the e±track 
channel. 

• Jet identification 

This uncertainty corresponds to the correction to 
simulated jet identification and quality requirement 
efficiencies described in Sec. IV El 

• Jet energy calibration 

This uncertainty includes contributions estimated 
for the jet energy scale and resolution corrections 
described in Sec. IV El 

• Trigger simulation 

Uncertainties on the fits to the trigger efficiencies 
for each object discussed in Sec. IIVI were propa- 
gated to estimate event triggering systematics. 

• Background estimation 

For background estimates using Monte Carlo simu- 
lations, normalization uncertainties were calculated 
using theoretical and/or experimental uncertainties 
in the corresponding product of production cross 
section and decay branching ratios. For instance, a 
systematic uncertainty of 35% is associated to the 
normalization of diboson+jets background, taken 
very conservatively as the difference between LO 
and NLO cross sections. For backgrounds esti- 
mated from data, systematic uncertainties have 
their ultimate origin in the statistical limitations 
of the relevant control samples. 

• tt tagging probability 

For the £+track channels, additional uncertainties 
associated with the b quark tagging probability in 
tt decays (Sec. IV Cp were included. The dominant 
sources of uncertainty arise from the method used 
to extract the b tagging efficiency in data and from 
the limited statistics in the heavy flavor enriched 
data samples. 

• Background tagging probability 

For the £±track channels, uncertainties in the 
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b tagging probabilities for the background pro- 
cesses (Sec. IVIEp were also taken into account. 
These originated from limited statistics in the 
background-enriched data samples and observed 
dependence of the event tagging probability on 
For the IF+jets events, where the tt contamination 
was subtracted assuming a cross section of 7 pb, 
the effect of varying the tt cross section between 5 
and 9 pb was propagated to the final result. 

• Luminosity 

The integrated luminosity corresponding to each of 
the data samples used by the analyses has a frac- 
tional uncertainty of 6.1% [HI]. 

• Other uncertainties 

Statistical uncertainties related to the sizes of 
Monte Carlo and data samples used independently 
for each channel are uncorrelated between them. 
These and other less important uncertainties are 
combined in this category. 

For each channel, an uncertainty on the cross section 
was obtained for each independent source of systematic 
uncertainty by varying the source by ±lcr and propagat- 
ing the variation into both background estimates and the 
signal efficiency. A new likelihood function was derived 
for each such variation to give a new optimal cross sec- 
tion. The resulting variations in the central value of the 
cross section are presented in Table IXVll 



C. Combined Cross Section 

Calculation of the combined estimate of the tt pro- 
duction cross section using all of the results presented in 
Table IXVI is complicated by the fact that some of the 
selection criteria are correlated. Specifically, the ee crite- 
ria are correlated with those of the e+track analysis, and 
the /i/i and /i+track criteria overlap. To account for this, 
we apply a BLUE technique (i.e., Best Linear Unbiased 
Estimate) [27[ . 

The correlations between the top quark pair selection 
efficiencies of the non-orthogonal analysis pairs were esti- 
mated with psuedo-experiments drawn from Monte Carlo 
samples. The ee-e+track and /z/x-/x+track correlations 
were found to be 43% and 47%, respectively. The use 
of b tagging in the £+track selections resulted in negligi- 
ble correlations between the backgrounds surviving each 
channel's selection. 

Correlations between the systematic uncertainties of 
each of the analyses were also included in the combina- 
tion. These were taken as 100%, —100%, or 0%, as appro- 
priate. Furthermore, all asymmetric uncertainties were 
made symmetric by averaging their positive and negative 
values. The combined cross section was derived using an 
iterative process. The combination in each iteration step 
was performed using the expected statistical and system- 
atic uncertainties evaluated at the cross section obtained 
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FIG. 9: The dependence of the measured cross section on the 
top quark mass compared to the theoretical prediction 0, [1[ . 



in the previous iteration step. The use of expected un- 
certainties avoids over-weighting the results of downward 
fluctuations. It was verified that the result of the itera- 
tive process was independent of the cross section input 
to the first iteration. The calculation was repeated until 
the result was stable to within 0.01% between iterations. 
The resulting combined cross section is 



o-ti 



= 7.4 ± 1.4 (stat) ± 0.9 (syst) ± 0.5 (lumi) pb (20) 

for a top quark mass of 175 GeV. There were a total of 
four degrees of freedom in the combination, and the x 2 
of the result is 1.6. Table IXVIII lists the relative weight of 
each analysis channel's result in the combination. The e/z 
measurement dominates the result, with the two £+track 
results entering with the next most significant weights. 
Table IXVIIII presents the contribution of each individ- 
ual systematic uncertainty described in Sec. lVIIBl to the 
total error on the result. 

The dependence of the result on the top quark mass 
was computed using parameterizations of each channel's 
selection efficiency as a function of mt- For a set of as- 
sumed masses, the individual channel results and their 
combination were recalculated using the appropriate ef- 
ficiency. The result is shown in Fig. [5] For values of m t 
between 170 GeV and 180 GeV, the value of the mea- 
sured cross section as a function of top quark mass is 
approximated by: 

a t t(m t ) = 7.4pb-0.1^ x {m t - 175 GeV). (21) 

For the current Tevatron average of top quark mass of 
170.9 GeV [28|, the resulting value of the cross section is 
7.8 ±1.8(stat+syst) pb. 
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TABLE XVI: A summary of the effects of individual systematic uncertainties on each channel's measured cross section. Quan- 
tities are presented in percent change from the central values presented in Table IXVl 



Source 


e/i 


ee 




e+track 


/K+track 


PV identification 


0.5 


—0.5 


0.5 


—0.5 


0.7 


—0.5 


0.6 


—0.6 


0.7 


—0.7 


Lepton identification 


6.3 


-5.9 


9.1 


-8.2 


3.8 


-3.5 


4.3 


-4.2 


7.8 


-7.2 


Track reconstruction 


4.2 


-4.6 


0.1 


-0.1 


8.8 


-10.3 


3.4 


-3.3 


4.5 


-4.0 


Jet identification 


4.8 


-3.5 


8.0 


-4.0 


13.2 


-8.6 


2.8 


-6.3 


6.6 


-0.8 


Jet energy calibration 


7.1 


-6.1 


8.2 


-5.4 


22.0 


-21.7 


8.5 


-11.2 


10.0 


-11.4 


Trigger 


10.2 


-5.7 


7.5 


-1.9 


5.5 


-4.0 


2.9 


-2.3 


5.5 


-4.4 


Bkg estimation 


4.7 


-3.7 


2.1 


-2.2 


5.6 


-5.4 


2.5 


-2.5 


3.8 


-3.9 


tt tagging 


0.0 


-0.0 


0.0 


-0.0 


0.0 


-0.0 


4.0 


-3.8 


4.0 


-3.8 


Bkg tagging 


0.0 


-0.0 


0.0 


-0.0 


0.0 


-0.0 


7.0 


-7.0 


11.2 


-11.2 


Other 


2.4 


-2.3 


2.3 


-2.2 


2.6 


-2.4 


2.2 


-2.3 


2.6 


-2.7 


Total 


16.2 


-12.5 


16.7 


-11.2 


28.6 


-26.7 


13.9 


-16.5 


20.5 


-19.6 



TABLE XVII: Relative weight of each measurement in the 
combined cross section calculation. 



Channel 


Weight (%) 


e/i 


53 


e+track 


22 


p+track 


17 


ee 


4 




4 


All 


100 


TABLE XVIII: A summary of the effects of individual system- 


atic uncertainties on the combined cross section measurement. 


Source 


Uncertainty (pb) 


PV identification 


0.07 


Lepton identification 


0.41 


Track reconstruction 


0.09 


Jet identification 


0.30 


Jet energy calibration 


0.60 


Trigger 


0.39 


Bkg estimation 


0.22 


tt tagging 


0.11 


Bkg tagging 


0.19 


Other 


0.10 


Total 


0.94 



VIII. CONCLUSION 

We have measured the tt production cross section in pp 
collisions at = 1.96 TeV utilizing dilepton signatures 
in approximately 425 pb -1 of data collected with the DO 
detector. The result, for m t = 175 GeV, is 

a ti = 7.4 ± 1.4 (stat) ± 1.0 (syst) pb. (22) 



This is in good agreement with the theoretical prediction 
of 6.7tQg pb from the full NLO matrix elements and 
the re-summation of the leading and next-to-leading soft 
logarithms 0, Q . For the current Tevatron average of m t 
= 170.9 GeV, the corresponding value of the measured 
cross section is 7.8 ±1.8 (stat+syst) pb. 
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APPENDIX A: TRIGGER REQUIREMENTS 

Tables IXIXI IXXIlTl list the trigger conditions used for 
each analysis channel. A description of the triggering sys- 
tem and details regarding particle reconstruction in each 
trigger subsystem are available in Sec. HVl As beam con- 
ditions changed and delivered instantaneous luminosity 
increased, trigger conditions were changed to maintain 
event selection rates within operational limits. The ta- 
bles group triggers used together at the same time and 
present the total luminosity exposed to each grouping. 
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TABLE XIX: Trigger requirements used to collect data for the e/x analysis. Total integrated luminosity exposed to each trigger 
set is given in the last column. 



Level 1 


Level 2 


Level 3 


Integrated 


conditions 


conditions 


conditions 


luminosity (pb _1 ) 


1 e, Et > 5 GeV and 1 loose fj, 


none 


1 loose e, E T > 10 GeV 


130.2 


1 e, Et > 6 GeV and 1 loose p, 


none 


1 loose e, E T > 12 GeV 


243.8 






1 loose e, E T > 12 GeV 








OR 




1 e, Et > 6 GeV and 1 loose n 


1/x 


1 tight e,E T >7 GeV 
OR 

1 tight, track-matched e, Et > 5 GeV 


53.2 



TABLE XX: Trigger requirements used to collect data for the ee analysis. Total integrated luminosity exposed to each trigger 
set is given in the last column. 



Level 1 


Level 2 


Level 3 


Integrated 


conditions 


conditions 


conditions 


luminosity (pb -1 ) 


2 e, E T > 10 GeV 


none 


1 loose e, St > 10 GeV 


23.3 


2 e, Et > 10 GeV 


none 


1 loose e, E T > 20 GeV 


120.2 


1 e, E T > 11 GeV 




2 loose e, E T > 20 GeV 




OR 




OR 




2 e, E T > 6 GeV 


none 


1 loose and 1 tight e, _Et > 15 GeV 


252.2 


OR 








2 e, 4 > 9 GeV, £ T > 3 GeV 








1 e, S T > 11 GeV 




2 loose e, S T > 20 GeV 




OR 




OR 




2 e, E T > 6 GeV 


2 e, ££t > 18 Ge y 


2 tight e, £ T > 8 GeV " 


49.8 


OR 




OR 




2 e, £ T > 9 GeV, E T > 3 GeV 




1 loose and 1 tight e, Et > 15 GeV 





"For part of the data, a 10 GeV Et requirement was used. 
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TABLE XXI: Trigger requirements used to collect data for the \xjx analysis. Each triggering regime has both a single muon and 
a dimuon requirement, each of which ties together trigger conditions at all three levels. Total integrated luminosity exposed to 
each trigger set is given in the last column. 



Multiplicity Level 1 

conditions 


Level 2 
conditions 


Level 3 
conditions 


Integrated 
luminosity (pb _1 ) 


1 1 tight /i 


1 M, 


> 3 GeV 1 track, p T > 10 GeV 


59.5 


2 2 loose /i 


1 


A* 


none 


1 1 tight /i 


1 M, PT 


> 3 GeV 1 track, p T > 10 GeV 




2 loose /i 

2 


1 


A 4 


1 /tt, p T >15 GeV 
OR 

1 track, p T > 10 GeV 


66.5 


1 1 tight /i 


1 A*, PT 


> 3 GeV 1 track, p T > 10 GeV 




2 loose /i 

2 


1 


A 4 


1 A*, Pt > 6 GeV 
OR 

1 track, pr > 5 GeV 


243.8 


1 1 tight /i and 1 track, p T > 10 GeV 1 At, pr 


> 3 GeV 1 Pt > 15 GeV 




2 loose (i 

2 


1 


A* 


1 A*, Pt > 6 GeV 
OR 

1 track, pt > 5 GeV 


51.5 


TABLE XXII: Trigger requirements used to collect data for the e+track analysis. For some periods of data collection, a logical 
OR of multiple requirements was used. Each requirement tied together trigger conditions at all three levels. Total integrated 
luminosity exposed to each trigger set is given in the last column. 


Level 1 
conditions 


Level 2 
conditions 




Level 3 
conditions 


Integrated 
luminosity (pb _1 ) 


1 e, E T > 10 GeV 

AND 

2 jets, E T > 5 GeV 


1 e, £ T > 10 GeV 
AND 

2 jets, £ T > 10 GeV 




1 tight e, E T > 15 GeV 
AND 

2 jets, E T > 15 GeV 


127.8 


1 e, E T > 11 GeV 


none 




1 tight e, E T > 15 GeV 
AND 
2 jets, £t > 20 GeV 




1 e, E T > 11 GeV 

OR 

2 e, E T > 6 GeV 

OR 

2 e, E T > 3 GeV, 1 e, E T > 9 GeV 


none 




1 tight e, E T > 20 GeV 


244.0 


1 e, E T > 11 GeV 

OR 

2 e, E T > 6 GeV 

OR 

2 e, E T > 3 GeV, 1 e, E T > 9 GeV 


none 




1 loose e, £t > 50 GeV 




le,£ T > 11 GeV 


lfi.Br> 15 GeV 




1 tight e, £t > 15 GeV 
AND 
2 jets, E T > 20 GeV 

AND 
1 jet, E T > 25 GeV 




1 e, E T > 11 GeV 

OR 

2 e, S T > 6 GeV 

OR 

2 e, £ T > 3 GeV, 1 e, E T > 9 GeV 


1 e, E T > 15 GeV 




1 tight e, E T > 20 GeV 


53.7 


1 e, E T > 11 GeV 


1 e, E T > 15 GeV 




1 loose e, £t > 50 GeV 





OR 

2 e, E T > 6 GeV 
OR 

2 e, E T > 3 GeV, 1 e, g T > 9 GeV 
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TABLE XXIII: Trigger requirements used to collect data for the /i+track analysis. For some periods of data collection, a logical 
OR of multiple requirements was used. Each requirement tied together trigger conditions at all three levels. Total integrated 
luminosity exposed to each trigger set is given in the last column. 



Level 1 


Level 2 


Level 3 


Integrated 


conditions 


conditions 


conditions 


luminosity (pb _1 ) 


1 loose fj, 


1 H 


1 jet, £ T > 20 GeV 




AND 






131.5 


1 jet E T > 5 GeV 








1 loose fi 
AND 


1 A 4 
AND 


1 jet, E T > 25 GeV 


244.0 


1 jet, E T > 3 GeV 


1 jet, E T > 10 GeV 




1 tight fi 


1 n, p T > 3 GeV 


1 track, p T > 10 GeV 




1 tight fi 


1 M 


1 jet, E T > 25 GeV 




AND 


AND 






1 jet, E T > 5 GeV 


1 jet, E T > 8 GeV 






1 loose n w. trackmatch 


none 


1 track, p T > 10 GeV 




AND 




OR 




1 track, p T > 10 GeV 




1 n, p T > 15 GeV 


30.3 


1 tight fi 
AND 


1 fi, pt > 3 GeV 


1 pr > 15 GeV 


1 track, p T > 10 GeV 








1 loose /i 
AND 


1 A* 


1 /i, pt > 15 GeV 




1 track, p T > 10 GeV 








1 tight fi 
AND 


1 M 
AND 


1 /tt, p T > 3 GeV 
AND 




1 jet, E T > 5 GeV 


1 jet, E T > 8 GeV 


1 jet, E T > 25 GeV 




1 loose n w. trackmatch 


none 


1 track, p T > 10 GeV 




AND 




OR 




1 track, p T > 10 GeV 




1 At, Pt > 15 GeV 


16.0 


1 tight fi 
AND 


1 Pt > 3 GeV 


1 At, Pt > 15 GeV 


1 track, p T > 10 GeV 








1 loose /i 


1 A* 


1 At, Pt > 15 GeV 





AND 

1 track, p T > 10 GeV 



